
Fig. 1 (Prior Art) 



Step 10. Load the current CB (cone beam) projection into computer 
memory. Suppose that the mid point of the CB projections currently 
stored in memory is yis^) . 



Steps 20 and 30, Depending on where the x-ray source is located (on 
L or onC), choose an appropriate family of filtering lines. Then, choose 
a discrete set of lines from the family. 



Step 40. Preparation for filtering. Parameterize points on the said 
filtering lines selected in Steps 20 and 30 by polar angle ;^ . Compute the 
derivative of the CB data (d / dq)D^(yiq), fi) l^^^at points /Son the said 
lines that correspond to discrete values of the polar angle. 
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Step 50. Filtering For each filtering line convolve the data for that line 
computed in Step 40 with filter 1 / sin ;^ . 
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Step 60. Back-projection. For each reconstruction point x , back- 
project the filtered data found in Step 50 according to equation (10). 
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Step 70. Go to Step 10, unless there 
process or image reconstruction at a 
completed. 


are no new CB projections to 
11 the required points x have been 



Reconstructed 
image 



Fig. 7 
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Step 20. Finding families of lines for filtering. 
It is assvuned the x-ray source is located on the line L . 



StcD 21. Choose a discrete set of values of the parameter s, inside 
the interval [^r",.?,"""] . 
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Steo 22. For each chosen in Step 21 find a line tangent to the 
projected circle C. 
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StcD 23. The collection of lines c 
required set of lines. 


onstructed in Step 22 is the 



step 30 
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Step 40. Preparation for filtering 

Step 41. If the x-ray source is located on the line Z , fix a fihering 
line Ljj^ e_ , . If the x-ray source is on the circle C , fix Ljj^ g_ 2 • 



Step 42. Parameterize points on the said line by polar angle / in 
the plane through yisQ) and Z,^, . 

Step 43, Choose a discrete set of equidistant values /j that will 
be used later for discrete filtering in Step 50. 



Step 44. For each yj find the unit vector fij which points fi-om 
y^s^) towards the point on I^, that corresponds to /j. 



Step 45, Using the CB projection data D^(y(qX®) for a few 
values of q close to Sq find numerically the derivative 
(d/dq)D^(y(q%e) |^^,^ for all 0 = /?, . 



Step 46, Store the computed values of the derivative in computer 
memory. 



Step 47. Repeat Steps 41-46 for all lines Ljj, identified in Step 20. 
This way we will create the processed CB data ^(s^./Jj) 
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Step 50 



Step 50. Filtering 

Step 51. Fix a filtering line L^^ . If the x-ray source is located on 
the line L , take Lj^^ . If the x-ray source is located on the 
circle C , take I^, €_ ^ • 



Step 52. Compute FFT of the values of the said processed CB data 
computed in Step 40 along the said line. 



Step 53. Compute FFT of the filter l/sin/ 



Step 54. Multiply FFT of the filter 1/sin;^ (the resuU of Step 53) 
and FFT of the values of the said processed CB data (the resuh of 
Step 52). 



Step 55. Take the inverse FFT of the resuh of Step 54. 



Step 56. Store the result of Step 55 in computer memory. 
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Step 57. Repeat Steps 51-56 for all lines in the said family of 
nines. This will give the filtered CB data <D(5o,y5,) . 



Step 60 



step 60. Back-projection 



Step 61 . Fix a reconstruction point x , which represents a point inside 
the patient where it is required to reconstruct the image. 

▼ 

Step 62. If Sq € /,(x)u/2(x), then the said filtered CB data affects the 

image at x and one performs Steps 63-68. If «s /, (x) u /jC^c) , then the 

said filtered CB data is not used for image reconstruction at x . In this 
case go back to Step 61 and choose another reconstruction point. 

Step 63. Find the projection x of x onto the detector plane DP^Sq) 
and the unit vector /9(5'o,jc), which points from yiso) towards x. 

▼ 

Step 64. Identify filtering lines L^, g_, orZ^, g_2 (depending on 

where the x-ray source is located) and points on the said lines that are 
close to X . This will give values of O(^o,/?^) for /3j close to J3(Sq,x) . 

Step 65. With interpolation estimate the value of ^{Sq,/3(Sq,x)) from 
the said values of <!>(Sq,/3j) for J3j close to /3{Sq,x), 

Step 66. Compute the contribution from the said filtered CB data to the 

image being reconstructed at the point x by multiplying 0(^0, /?(5o,x)) 

by -Sf^(Sf^,x)/(27r^ \x-y{sQ)\) (see (8) for definition of ^^^(^o*^))- 
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Step 67. Add the said contribution to the image being reconstructed at 

the point x according to a pre-selected scheme for approximate 

evaluation of the integral in equation (10). 

▼ 

Step 68. Go to Step 61 and choose a different reconstruction point x . 



